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Abstract 

We investigate the possibility of electroweak phase transition in the minimal 
supersymmetric standard model (MSSM) with an extra U(l)' . This model has two 
Higgs doublets and a singlet, in addition to a singlet exotic quark superfield. We 
find that at the one-loop level this model may accommodate the electroweak phase 
transitions that are strongly first-order in a reasonably large region of the parameter 
space. In the parameter region where the phase transitions take place, we observe 
that the lightest scalar Higgs boson has a smaller mass when the strength of the 
phase transition becomes weaker. Also, the other three heavier neutral Higgs bosons 
get more large masses when the strength of the phase transition becomes weaker. 



1 



I. INTRODUCTION 



The baryon asymmetry of the universe can be dynamically generated during the evolution 
of the universe, if the mechanism of baryogenesis satisfies the three Sakharov conditions 
[1]. The three Sakharov conditions are: the presence of baryon number violation, the 
violation of both C and CP, and a deviation from thermal equilibrium. It is known that 
the universe can escape out of the thermal equilibrium by means of electroweak phase 
transition, which should be strongly first-order in order to ensure sufficient deviation from 
thermal equilibrium to generate the baryon asymmetry that is observed today. However, 
it has been already recognized that the Standard Model (SM) has some difficulty to realize 
the desired electroweak phase transition. The present experimental lower bound on the 
mass of the SM Higgs boson does not allow the electroweak phase transition to be strongly 
first-order [2, 3]. The electroweak phase transition is weakly first-order or higher order in 
the SM. Thus, the SM is inadequate to generate sufficient baryon asymmetry. Moreover, 
the amount CP violation in the Cabibbo-Kobayashi-Maskawa (CKM) matrix is too small 
to account for the baryon asymmetry of the observed universe [4]. 

Consequently, new physical models beyond the SM have extensively been studied for 
the possibility of reasonable explanation of the baryon asymmetry of the universe. Espe- 
cially, the low energy supersymmetric models have been studied widely within the context 
of electroweak baryogenesis [5-7]. The simplest supersymmetric model that includes the 
SM is the minimal supersymmetric standard model (MSSM), which possesses in its su- 
perpotential the \i term that accounts for the mixing between two Higgs doublets. The 
yU parameter, which has the mass dimension, causes some problem with respect to its 
energy scale [8]. Several possibilities have been investigated in the literature to solve the 
so-called fi problem [9-12]. Introducing an additional U(l)' to the MSSM is one of the 
plausible explanations for the /i problem of the MSSM. 

The MSSM with an extra U(l)' can not only solve the \i problem but we will show 
that it can also overcome the difficulties that the SM encounters when the SM tries to 
satisfy the Sakharov conditions. This model can accommodate sufficient CP violation, 
because it possesses other sources of CP violation besides the CKM matrix. It is possible 
to realize the explicit CP violation in this model by means of complex CP phases arising 
from the soft SUSY breaking terms [12]. 

Then, it is the purpose of this paper to show that this model indeed allows the strongly 
first-order electroweak phase transitions such that it can successfully explain the baryo- 
genesis. The characteristics of the electroweak phase transitions are determined essen- 
tially by the temperature-dependent part of the Higgs potential. We construct the full 
temperature-dependent Higgs potential at the one-loop level, and examine if the elec- 
troweak phase transition may be strongly first-order. Two methods are employed for the 
construction of the temperature-dependent Higgs potential. One method assumes that 
the critical temperature at which the electroweak phase transition occurs is relatively 
high, thus the temperature-dependent effective potential is approximated by retaining 
only terms proportional to T 2 , whereas the other method carries out numerically exact 
integrations of the temperature-dependent effective potential. The thermal effects of par- 
ticles whose masses are comparatively smaller than the critical temperature are included 
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at the one-loop level in the former method, whereas the particle content is different in the 
latter method. 

Either way, we obtain almost the same physical results. Unlike the MSSM, this model 
allows a strongly first-order electroweak phase transition in a wide region of the parame- 
ter space, and the first-order electroweak phase transition can be strong enough without 
requiring a light stop quark. An interesting behavior of this model with respect to the 
strongly first-order electroweak phase transition is that the mass of the lightest neutral 
Higgs boson becomes larger when the phase transition gets stronger. On the other hand, 
the masses of the other three neutral Higgs bosons become smaller when the phase tran- 
sition gets stronger. 



II. ZERO TEMPERATURE 

The MSSM with an extra U(l)' accommodates in its Higgs sector two Higgs doublets 
Hi = (H®, Hi), H 2 = (H 2 , H®), and one Higgs singlet, S. In terms of these Higgs fields, 
the relevant part of the superpotential of this model may be written as 

W « h t QH 2 t c R + h b QHib c R + h k SD L D R - \SH^eH 2 , (1) 

where we take into account only the third generation: t c R and b c R are, respectively, the 
right-handed singlet top and bottom quark superfields, D R is the right-handed singlet 
exotic quark (a vector-like down quark) superfield, Q is the left-handed SU (2) doublet 
quark superfield of the third generation, and Dl is the left-handed singlet exotic quark 
superfield. Further, h t , h\, and hk are, respectively, the dimensionless Yukawa coupling 
coefficients of top, bottom, and exotic quark superfields, and e is an antisymmetric 2x2 
matrix with 6i 2 = 1. 

From the superpotential, at zero temperature, we can construct the Higgs potential 
at the tree level, which may be read as 

V = V F + V D + V s , (2) 

where 

V F = \\\ 2 [(\Hi\ 2 + \H 2 \ 2 )\S\ 2 + \H^eH 2 \ 2 }, 

V D = 9 4(H\*Hi + Hl*H 2 ) 2 + 9 4(\Hi\ 2 - \H 2 \ 2 ) 2 

+ 9 l(Q 1 \Hi\ 2 + Q 2 \H 2 \ 2 + Q 3 \S\ 2 ) 2 , 
V s = m\\Hi\ 2 + m\\H 2 \ 2 + ml\S\ 2 - [\A x (H?eH 2 )S + H.c] , (3) 

where a denotes the three Pauli matrices, gi, g 2 , and g[ are the U(l), 577(2), and U(l)' 
gauge coupling constants, respectively, Qi, Q 2 , and Q 3 are the U(l)' hypercharges of Hi, 
H 2 , and S, respectively, and mf (i = 1, 2, 3) are the soft SUSY breaking masses. In the 
Higgs potential, A and A\ may in general be complex numbers. However, they will be 
assumed to be real in the subsequent discussions, as we do not consider CP violation in 
the Higgs sector. The soft masses are also assumed to be real, without loss of generality, 
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and they are eventually eliminated by imposing minimum conditions with respect to the 
neutral Higgs fields, The gauge invariance of the superpotential under of U(l)' requires 
that the three U(l)' hypercharges should satisfy Qi + Q 2 + Q3 = 0. 

The above Higgs potential at the tree level would allow the three neutral Higgs fields 
H±, H®, and S to develop the vacuum expectation values (VEVs) i>i(0), 1*2(0), and s(0), 
respectively. Remark that these VEVs are obtained at zero temperature. However, for 
simplicity, we omit the temperature dependence of these VEVs until next section where 
we take into account the finite temperature effect. 

The tree-level Higgs potential should now be corrected by the radiative one-loop effects. 
In SUSY models, the radiative corrections due to the top and stop quarks contribute 
most dominantly to the tree-level Higgs sector. Besides, if tan/5 = 1)2/1)1 is very large, 
the radiative corrections due to the bottom and sbottom quarks should also be included 
since they become no longer negligible. Furthermore, the radiative corrections due to the 
exotic quark and squark may be important since the Yukawa coupling of the exotic quark 
to the singlet field S can be large at the electroweak scale [11]. Therefore, we take into 
account all the contributions from the top, bottom, exotic quark sector to the tree-level 
Higgs potential. 

The one-loop radiative corrections are evaluated by the effective potential method [13]. 
We assume that the squark masses are degenerate. Ignoring the mixings in the masses of 
the squarks [14], the one- loop effective potential is given by 



167T 2 
l=t,b,k 1U/I 



3 , (m 2 + MT 

2 + log {-mi 
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(4) 



where t, b, and k, respectively are top, bottom, and exotic quark fields including the 
corresponding squark fields, M. t = ht\H 2 \, M.b = h b \Hi\, M. k = h k \S\ are the field- 
dependent quark masses, and fh is the soft SUSY breaking mass, which is assumed that 
m = 1000 GeV > m q (q= t, b, or k). 

The Higgs sector of the present model consists of six physical Higgs bosons: a pair 
of charged Higgs boson, one neutral pseudoscalar Higgs boson, and three neutral scalar 
Higgs bosons. The tree-level mass of the charged Higgs boson is given by 

™ c ± = m w - A v + . , (5) 

sin _ lj 



where v = ywf + v\ = 175 GeV and = g\v 2 /2 is the squared mass of the W boson. 
At the tree level, the mass of the charged Higgs boson might be either smaller or larger 
than the W boson mass. 

The tree-level mass of the neutral pseudoscalar Higgs boson is given by 

sin 2a 

where tan a = {v/2s) sin 2/3 implies the splitting between the electroweak symmetry break- 
ing scale and the extra U(l)' symmetry breaking scale. Note that these tree-level masses 
of both the neutral pseudoscalar and the charged Higgs bosons do not receive any radiative 
corrections, because the squark masses are degenerate. 
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The tree-level squared masses of the three neutral scalar Higgs bosons are considerably 
affected by the radiative corrections. Their squared masses at the one-loop level are given 
as the eigenvalues of the 3x3 one-loop level mass matrix, whose elements may be written 

as 



M u 
M 22 
M 33 
M 12 
M 13 
M 23 

where m 2 z 
is defined as 



m z cos f3 + 2g 1 Q 1 v cos f3 + m A sm (3 cos a + f a ( m b) 
m 2 z sin 2 (3 + 2g 2 Q\v 2 sin 2 (3 + m 2 A cos 2 (3 cos 2 a + f a { m Ti > 

o '2/S2 2 1 2-2 1 / 2 \ 

2g x Q z s +m A sm a + f a (m k ), 

91Q1Q2V 2 sin 2/3 + (A 2 f 2 — m 2 z /2) sin 2/3 — cos /3 sin /? cos 2 a , 
2g 2 QiQ 3 vs cos /3 + 2\ 2 vs cos /3 — sin (3 cos a sin a , 
2giQ2Q?,vs sin /3 + 2A 2 us sin (3 — m 2 A cos (3 cos a sin a , 



(7) 



+ 5 , 2) t ' 2 /2 is the squared mass of the Z boson, and the function f a (jn 2 q ) 



' ~ 2 
m 



+ 



m 2 + m 2 



(m 2 + m 2 ) 2 



We assume that the masses of three scalar Higgs bosons Si are sorted such that m Sl < 



III. FINITE TEMPERATURE 



Now, let us study the temperature dependence of the Higgs potential in order to inves- 
tigate the nature of the electroweak phase transition in the MSSM with an extra U(l)'. 
We evaluate Vp, the temperature-dependent part of the Higgs potential at the one-loop 
level, using the effective potential method. It is given as [15] 



V r 



v UlTi 

^ 2tt 2 

l=B,F Z " 



dx x 2 log 



1 ± exp (^- v /x 2 + m 2 (0 i )/T 2 ) 



(9) 



where B and F stand for bosons (t, b, and k) and fermions (t, b, and k), and n t = nt, = 
Uk = —12 and ni = n~ h = nr = 12. The negative sign is for bosons and the positive sign 
is for fermions. Thus, the full Higgs potential at finite temperature at the one-loop level 
is given by 



V(T) = V + Vx + V T 



(10) 



For numerical analysis, we need to set the values of the relevant parameters of the 
model. As in the previous section, the soft SUSY breaking mass is set as rh = 1000 GeV. 
The quark masses ar e set as m t = 175 GeV, = 4 GeV, and = 400 GeV. From 

t, b, k) yield the squark masses as mi 



these values, 



rrig 



m 2 



(q 



1015 GeV, 



= 1000 GeV, and mf. 



1077 GeV. 

Some caution should be taken for setting the values of Qt (i=l, 2, 3), the U{1)' 
hypercharges of the Higgs doublets and the Higgs singlet. In the MSSM with an extra 
£7(1)', the extra neutral gauge boson mass (mz') and the mixing angle (oczz') between the 
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two neutral gauge bosons (Z, Z') may impose strong constraints on the parameter values. 
For our numerical analysis, mz 1 is estimated to be larger than 600 GeV, and azz 1 smaller 
than 2 x 10~ 3 , for tan/5 = 3 and s(T = 0) = 500 GeV. Besides, as recent research has 
suggested [10], we impose the constraint of Q1Q2 > 0. Further, the U(l)' gauge invariance 
condition requires that Q3 = —(Qi + Qz)- 

In this paper, we define new charges Qi = g[Qi since Qi appear always together with 
g[. Then, one may establish the allowed area in the (Qi, (^-plane by imposing the above 
constraints. For tan/3 = 3 and s(T = 0) = 500 GeV, the result is shown in Fig. 1, where 
the small area near the point (Qi, Q 2 ) = (-1, 0) and the upper right corner of Fig. 1 
are the allowed areas. The hatched region is the excluded area. There are two specific 
points in Fig. 1, marked by a star (*) and a cross (+). The values of Qi and Q2 at 
the star-marked point correspond to the i/-model of E 6 gauge group realizations [11]. We 
would take the values of Q\ and Q2 at the cross-marked point, namely, (Qi,Q2) = (-1; 
-0.1), and hence Q3 =1.1. 

With these parameter values at hand, we would investigate the possibility of the 
strongly first-order electroweak phase transition by using two different ways. The first 
method is to retain only the dominant T 2 -proportional part from the high-temperature 
approximation of Vr, and to take account only those particles whose masses are relatively 
small [6]. The second method is to perform the integration in Vr in numerically exact way, 
and to consider only the contributions of top, bottom, and exotic quarks and squarks. 



1. Method A 

Let us start with the high temperature approximation of Vr, which is expressed as [3] 
_ v \ T'rnH<P.) , mKfr) ( rnHM Y 

+ U i 
i=t,b,k 

where log cp = 2.64 and logcg = 5.41. It is known that in the SM the high temperature 
approximation is consistent with the exact integration of Vp within 5 % at temperature 
T for mp/T < 1.6 and rns/T < 2.2, where mp and mg are respectively the fermion mass 
and the boson mass that participate in the potential. 

We select those terms that are proportional to T 2 in the above expression, which 
become most dominant at high temperature. Thus, we assume that the temperature at 
which the electroweak phase transition takes place is sufficiently high. We also assume 
that the U(l) and SU(2) gaugino masses Mi and M 2 in the chargino and neutralino 
sectors are very much larger than the other mass parameters. We take into account the 
thermal effects due to the Higgs bosons, W, Z, and the extra U(l) gauge boson in the 
boson sector, and t, b, k quarks, the lighter chargino, and the three light neutralinos in the 
fermion sector, because their masses are relatively small as compared with temperature, 
similarly to the analyses of previous articles [6]. Explicitly, the T 2 terms in the high 



T 2 m 2 (0,) Tmf(0 t ) mf(&) / m 2 (0 t ) \ 
24 12tt 64vr 2 ° g 1 c B T 2 ) 
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temperature approximation of Vp can be expressed as 
T 2 



21 



4m 2 + 4m 2 , 



2m 2 3 + (2g{ + 6g 2 2 + 6A + \H 2 



12A 2 |Sf 



+ \2 gi \Q{\HiV + Qtm 1 + Ql\SV) + 2g 2 Q x Q 2 (\E x 
+ 2g' 2 Q 2 Q 3 (\H 2 \ 2 + \S\ 2 ) + 2g' 2 Q 1 Q 3 (\H 1 \ 2 + \S\ 2 ) 
+ %; 2 (Qi + Q^iQ^l 2 + Q 2 \H 2 \ 2 + Q 3 \S\ 2 ) 
+Q(h 2 t \H 2 \ 2 + h 2 b \H l \ 2 + h 2 k \S\< 



(12) 



Now, the neutral scalar Higgs fields develop the temperature-dependent VEVs, t>i(T), 
v 2 (T), and s(T), which we will simply denote v\, v 2 , and s, respectively. In terms of 
these temperature-dependent VEVs, the vacuum at finite temperature is defined as the 
minimum of V(T) as 

(V(v u v 2 , s, T)> = (V ) + (V x ) + (V T ) , (13) 

where 

c 2 j_ „2 



(Vo) 

(Vi) 
(Vt) 



,2..2 , „2..2 , ^2„2 , ^1 t^ (u 2_ u 2 ) 2 + A 2 ( ^2 + t; 2 s 2 + t; 2 s 2 ) 



'2 

hijn 2 ) + / 6 (m 2 ) + f b (m 2 k ) , 

T 2 r 

— Ami + 4m 2 + 2m s + ( 2 9i + 6 ^2 + 6A 2 )(w 2 + u 2 ) + 12AV 



+ 12^ 2 (g 2 t; 2 + g 2 W 2 



2 ^Qls z ) + 2g 1 2 Q 1 Q 2 (vt + v 2 2 ] 



>V) + 2#j 2 ( 
+ 2g' 2 Q 2 Q 3 (v 2 2 + s 2 ) + 2g 2 Q l Q 3 (v 2 l + s 2 ) 
+8g' 2 (Q 1 + Q 2 )(Qivl + Q 2 v 2 2 + Q 3 s 2 ) + 6(h 2 t v 2 2 + h 2 h v\ + k 2 s 2 ] 
In the above expressions, the function /& is defined as 



3m^ 

16^ 



"3 , 


/m 2 + m 2 \~ 


2 + bg 









(14) 
(15) 



and the soft SUSY breaking masses at the one-loop level are given as 



ffl 



m. 



m 



cos 2(5 - A 2 (s(0) 2 + w(0) 2 sin 2 (3) + AA A s(0) tan (3 



m 



_9 1 2 Qi(Qiv(0) 2 cos 2 + Q 2 v(0) 2 sin 2 + Q 3 s(0) 2 ) - / c (m 2 (0)) 

cos 2/? - A 2 (s(0) 2 + v(0) 2 cos 2 0) + AA A s(0) cot /5 
' ^ 2 Q 2 (Qi^(0) 2 cos 2 (3 + Q 2 v(0) 2 sin 2 /? + Q 3 s(0) 2 ) - / c (m 2 (0)) 
A 2 ^(0) 2 + ^yt;(0) 2 A A sin2/3 



- ^ 2 Q 3 (0i^(0) 2 cos 2 (3 + g 2 t'(0) 2 sin 2 (3 + Q 3 s(0) 2 ) - / c (m 2 (0)) , (16) 
where fi(0), v 2 (ti), and s(0) are the VEVs evaluated at zero temperature in the preceding 
section, tan/? = t> 2 (0)/t>i(0), v(0) = y / t'i(0) 2 + t>2(0) 2 = 175 GeV, and the function f c is 
defined as 

,2 



3/z 2 m 2 



'rh 2 + mV 



m„ 



Now, let us determine the critical temperature at which the electroweak phase tran- 
sition takes place. In our analysis, the critical temperature is defined by a temperature 
at which (V(T)) has two distinct minima with equal value, that is, a pair of degenerate 
vacua. In order to have a pair of degenerate vacua, the potential (V(T)) should satisfy 
the minimum condition of 

= 2m\s -2\A x v l v 2 + 2\ 2 (vl + v 2 2 )s 

+ 2g' 2 Q 3 s(Q 1 v 2 1 + Q 2 v\ + Q 3 s 2 ) + 2h 2 k m k f c (m 2 k ) 

+ ^s[2AX 2 + 2Ag' 2 Q\ + 20g' 2 Q 3 (Q 1 + Q 2 ) + 12k 2 } , (18) 

which is obtained by calculating the first derivative of the full effective potential at the 
finite temperature with respect to s. 

For given parameter values at given temperature, one may solve the above minimum 
condition to express s in terms of the other two VEVs, V\ and v 2 - Then, by substituting 
s into (V(vi,V2,s,T)), one may obtain (V(vx,V2,T)) which depends only on v± and v%. 
By inspecting the shape of (V(v±, V2, T)) on the (vi, ^-plane for given parameter values 
at given temperature, we may determine whether it possess a pair of degenerate vacua or 
not. 

In Fig. 2, the equipotential contours of (V(vi,V2,T)) are plotted on the (i>i,i>2)- 
plane, where the parameter values are set as tan/5 = 3, A = 0.8, s(0) = 500 GeV, 
rriA = 1830 GeV, and the temperature is set as T = 100 GeV, which is actually the 
critical temperature T c . One can easily spot two distinct minima of (V(v\,V2,T)) on 
the (v i, w 2 )-plane, namely, one at (0,0) and the other at (275,640) GeV. The phase of 
the state is symmetric at the minimum point (0,0) on the (fi, t^)-plane, whereas it is 
broken at (275, 640) GeV. The electroweak phase transition may take place from (0, 0) to 
(275,640) GeV on the (v±, t>2)-plane, which is evidently discontinuous and therefore it is 
first-order. 

The distance on the (vi, u 2 )-plane between the two minima of (V(vi,v 2 ,T)), defined 
as v c , determines the strength of the electroweak phase transition. The electroweak phase 
transition is said to be strong if v c /T c > 1, and weak otherwise. In Fig. 2, the distance is 
calculated to be 

v c = ^(275 - 0) 2 + (640 - 0) 2 = 696 (GeV) . (19) 

In Fig. 2, the strength of the electroweak phase transition is about v c /T c = 6.9, which 
definitely tells that the electroweak phase transition is a strong one. Therefore, the 
particular parameter values set for Fig. 2 yields an electroweak phase transition which 
is first-order as well as strong. Note that v c does not depend on s, that is, we need not 
to know the values of s at the two minima to calculate v c . Actually, v c is the VEV at 
the broken phase. The masses of the neutral scalar Higgs bosons at zero temperature 
for the parameter values of Fig. 2 are obtained as ms 1 = 56 GeV, ms 2 = 807 GeV, and 
m Ss = 1827 GeV. 

We repeat the above job of analysis, varying the values of the relevant parameters. We 
find that there are a large number of sets of parameter values that allow strongly first-order 
electroweak phase transitions. Thus, the MSSM with an extra U(l)' may accommodate 
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TABLE 1: Some sets of A and that allow strongly first-order electroweak phase 
transitions in the MSSM with an extra U(l)', obtained by Method A. The values of other 
parameters are fixed as tan/3 = 3, s(0) = 500 GeV, m = 1000 GeV, and T c = 100 
GeV. The pair of numbers in the third column are the coordinates of the broken-phase 
minimum of (V(vi, v 2 , T)). The coordinates of its symmetric-phase minimum is (0,0) 
for all sets. The three numbers in the fourth column are the masses of S±, S 2 , and £3, 
respectively. The number in the last column is the strength of the first-order electroweak 
phase transition. 



A 


m A (GeV) 


( Vl ,v 2 ) (GeV) 


m Sl , m S2 , 7715.3 (GeV) 


v c /T c 


0.1 


478 


(1750, 1650) 


120, 524, 792 


26 


0.2 


675 


(1400, 1500) 


118, 674, 796 


23 


0.3 


900 


(1200, 1400) 


112, 786, 908 


18 


0.4 


1109 


(870, 1200) 


104, 792, 1112 


15 


0.5 


1306 


(600, 1000) 


93, 796, 1307 


12 


0.6 


1486 


(430, 850) 


82, 800, 1485 


8 


0.7 


1660 


(340, 700) 


70, 803, 1658 


7 


0.8 


1830 


(275, 640) 


56, 807, 1827 


6.9 



the desired phase transitions for a wide region in its parameter space. Some of the results 
are listed in Table 1, where tan/3 = 3, s(0) = 500 GeV, and T = 100 GeV are fixed as 
the values set in Fig. 2, whereas A and have different values. The set of numbers in 
the last row of Table 1 is the numerical result of Fig. 2. 

Every set of numbers in each row of Table 1 gives (V(vi,v 2 ,T)) a pair of degenerate 
minima, the minimum of symmetric phase at (0,0) on the (v 1, ^2) -plane, and the one of 
broken phase at a different point on the (vi, t^-plane as given in Table 1. The electroweak 
phase transition is strongly first-order. One may easily observe in Table 1 that, as the 
value of A increases, a larger value of rriA allow desired phase transitions. On the other 
hand, the strength of the phase transition is reinforced if the value of A decreases. 

The masses of the neutral scalar Higgs bosons exhibit some interesting behavior. For a 
larger value of m^, both S 2 and S3 have also larger masses whereas Si has a smaller mass. 
The tendency is that the strength of the phase transition is reinforced if ms 1 increases 
and if m^, ms 2 , and ms 3 decrease. In the SM, the strength of the first order electroweak 
phase transition decreases if its single Higgs boson mass is increased. Also, in the MSSM, 
we have a weaker phase transition if the lighter one of its two scalar Higgs bosons has a 
larger mass. In this regard, the tendency of our model is opposite to those of the SM or 
the MSSM. One can see that this strange behavior also occurs in some parameter region 
of a non-minimal SUSY model, as shown in Fig. 3 of Ref. [7]. 



9 



2. Method B 



The second method evaluates Vt by exact integration to obtain the temperature-dependent 
full potential V(T) at one-loop level, where the thermal effects of top, bottom, and exotic 
quarks and squarks are taken into account. The thermal effects of the gauge bosons can 
be a help for strengthening the first-order electroweak phase transition, but we would 
omit them, since the strength of the phase transition is already strong enough. 

This method starts with the exact integral expression for (Vt) after replacing the 
neutral Higgs fields by their VEVs as 



(V T ) 



6i r°° 
- Y] —7- / dx x 2 log 



l=i,b,k 



7l z JO 



1 — exp 




mf(v u v 2 ,s) 



rp2 



1 + exp —\x 2 + 



m 2 + mf (v\, v 2 , s) 



(20) 



which is different from (Vt) of Method A, while (Vo) and (Vi) are the same as those of 
Method A. From the full (V(T)) = (Vo) + (Vi) + (Vt), we obtain a minimum condition 
for degenerate vacua as 

= 2mls-2\Axv 1 v 2 + 2\ 2 (v 2 + v 2 2 )s + 2g[ 2 Q 3 s(Qiv 2 + Q 2 vi + Q3S 2 ) 
+ ^h 2 k m k f c (m 2 k ) 

3T 2 r 00 » 2h 2 k s exp(- Jx 2 + m 2 k /T 2 ) 

H — / dx x 



IT 2 Jo 



x 2 + m 2 k /T 2 1 + exp(- Jx 2 + m 2 k /T 2 ' 



7T 2 JO 



dx x' 



X 



2h 2 k sexp{-Jx 2 + (m 2 + m 2 )/T 2 ) 



x 2 + (m 2 + ml) /T 2 1 + exp(-yx 2 + (m 2 + m 2 k ) /T 2 )) 



(21) 



where depends only on s and is independent from v\ and v 2 - 

Solving the above minimum condition is harder than solving the corresponding mini- 
mum condition of Method A. Nevertheless, we can solve it by using the bisection method 
to express s in terms of the other parameters. Then, eliminating s from (V(T)), we can 
obtain the expression for (V(vi,v 2 ,T)) which depends only on V\ and v 2 . Subsequent 
steps of numerical analysis are the same as the previous method. 

In Fig. 3, equipotential contours of (V(vi, v 2 , T)) obtained by the present method is 
plotted on the (vi,v 2 )-p\ane, where the parameter values are set slightly different from 
the previous method: tan (3 = 3, A = 0.8, s(0) = 500 GeV, m A = 1780 GeV, and T = 100 
GeV. The shape of the equipotential contours of Fig. 3 is almost the same as that of 
Fig. 2. One can see that there are two distinct minima in Fig. 3, just like Fig. 2: 
one at (0,0), and the other at (165,440) GeV on the (v\, u 2 )-plane, indicating that the 
phase transition is first order. The strength of the first-order phase transition is strong, 
since v c /T c = 4.7. The masses of the three scalar Higgs bosons are evaluated at zero 
temperature as ms 1 = 82 GeV, mj 2 = 804 GeV, and ms 3 = 1777 GeV. 
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TABLE 2: Some sets of A and that allow strongly first-order electroweak phase 
transitions in the MSSM with an extra U(l)', obtained by Method B. Other descriptions 
are the same as Table 1. 



A 


vrtA GeV 


{VlB,V2B) GeV 


m Si GeV 


Vc/J-c 


0.1 


462 


(1600, 1600) 


121, 468, 791 


22 


0.2 


663 


(1400, 1400) 


118, 662, 795 


19 


0.3 


885 


(1100, 1100) 


113, 785, 894 


15 


0.4 


1095 


(800, 1200) 


106, 792, 1098 


14 


0.5 


1287 


(680, 990) 


97, 796, 1288 


12 


0.6 


1457 


(400, 750) 


91, 799, 1456 


8 


0.7 


1620 


(300, 600) 


86, 801, 1618 


6 


0.8 


1780 


(165, 440) 


82, 804, 1777 


4.7 



Comparing Fig. 3 with Fig. 2, one may safely remark that Method A and Method 
B lead qualitatively the same results. Either method, whether (Vp) is calculated by 
direct integration or is simplified by high-temperature approximation, and whether the 
participating particles at the one-loop level are somewhat exhaustive or selective, we 
find that the MSSM with and extra U(l)' allows strongly first-order electroweak phase 
transitions for certain region in its parameter space. 

We repeat the numerical analysis by varying the parameter values, and some of the 
results are listed in Table 2. Like in Table 1, tan/3 = 3, s(0) = 500 GeV, and T = 100 
GeV are fixed, whereas A and wia are varied. The set of numbers in the last row of Table 
2 is the numerical result of Fig. 3. Comparing Table 2 with Table 1, one may notice that 
the numbers are slightly different from each other but the general behavior of the two 
tables is exactly the same. 

IV. DISCUSSIONS AND CONCLUSIONS 

We investigate the MSSM with an extra U(l)' if it could accommodate strongly first- 
order electroweak phase transitions to provide sufficient baryon asymmetry, for reasonable 
masses of scalar Higgs bosons. To do so, we need the temperature-dependent part of 
the Higgs potential at the one-loop level. Explicitly, its expression is obtained by two 
complementary methods: Method A employs high-temperature approximation and retains 
only the most dominant T 2 terms, and takes into account the thermal effects at the one- 
loop level of various participating particles. On the other hand, method B performs 
numerical integrations, and the thermal effects of top, bottom, and exotic quarks and 
squarks are accounted for. 

Both methods lead us to essentially the same conclusion: the strongly first-order 
electroweak phase transition is possible in the MSSM with an extra U(l)', for a wide 
region in its parameter space. The masses of the scalar Higgs bosons are obtained within 



11 



reasonably acceptable ranges. Accordingly, we may expect that the MSSM with an extra 
U(l)' can explain the baryon asymmetry of the universe. 

We remark that the MSSM with an extra U(l)' exhibits an interesting behavior with 
respect to the correlation between the strength of the phase transition and the Higgs 
boson masses. The MSSM with an extra U(l)' is opposite to the SM or to the MSSM 
in the sense that the mass of the lightest scalar Higgs boson increases when the strength 
of the strongly first-order electroweak phase transition becomes stronger. In the SM, its 
single Higgs boson has a larger mass when the strength of the first order electroweak 
phase transition decreases. In the MSSM, we also have a larger mass for the lighter one 
of its two scalar Higgs bosons when the phase transition becomes weaker. 
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FIGURE CAPTION 

FIG. 1. : The allowed area in the (Q\, Q 2 )-plane. For tan/3 = 3 and s(T = 0) = 500 
GeV, the small area near the point (Qi, Q2) = (-1, 0) and the upper right corner are the 
allowed areas, whereas the hatched region is the excluded area. There are two specific 
points, marked by a star (*) and a cross (+). The values of Q\ and Q 2 at the star-marked 
point correspond to the z/-model of E 6 gauge group realizations. The values of Q\ and Q 2 
at the cross-marked point are (Qi, Q 2 ) = (-1, -0.1), and hence Q3 =1.1. In our discussions, 
we choose this point. 

FIG. 2. : The plot of the equipotential contours of (V(vi,v 2 ,T)) on the (v±, u 2 )-plane, 
obtained by Method A. The parameter values are set as tan/3 = 3, A = 0.8, s(0) = 500 
GeV, rriA = 1830 GeV, and the temperature is set as T = 100 GeV, which is actually the 
critical temperature T c . Notice two distinct minima of (V(v\, v 2 , T)) on the (vi, i^-plane: 
(0, 0) where the phase of the state is symmetric, and (275, 640) GeV, where the phase 
of the state is broken. The electroweak phase transition may take place from (0, 0) to 
(275,640) GeV on the (v\, i>2)-plane, which is evidently discontinuous and therefore it is 
first order. The distance between the two minima is v c = 696 GeV, indicating that the 
strength of the first-order phase transition is strong (v c /T c > 1). The masses of the three 
scalar Higgs bosons are obtained as m Sl = 56 GeV, m S2 = 807 GeV, and ms 3 = 1827 
GeV. 

FIG. 3. : The plot of the equipotential contours of (V(vi,v 2 ,T)) on (vi, u 2 )-plane, ob- 
tained by Method B. The parameter values are set as tan/3 = 3, A = 0.8, s(0) = 500 
GeV, rriA = 1780 GeV, and T c = 100 GeV. The coordinates of two minima are: (0, 0) 
and (165,440) GeV. The distance between the two minima is v c = 470 GeV, thus the 
electroweak phase transition between the two minima is strongly first-order. The masses 
of the three scalar Higgs bosons are obtained as = 82 GeV, ms 2 = 804 GeV, and 
m S3 = 1777 GeV. 
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FIG. 1: The allowed area in the (Qi, Q2)-plane. For tan (3 = 3 and s(T = 0) = 500 
GeV, the small area near the point (Qi, Q2) = (-1, 0) and the upper right corner are the 
allowed areas, whereas the hatched region is the excluded area. There are two specific 
points, marked by a star (*) and a cross (+). The values of Q\ and Q2 at the star-marked 
point correspond to the z/-model of E$ gauge group realizations. The values of Q\ and Q2 
at the cross-marked point are (Qi, Q2) = (-1, -0.1), and hence Q3 =1.1. In our discussions, 
we choose this point. 
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FIG. 2: The plot of the equipotential contours of (V(vi,v 2 ,T)) on the (vi, i> 2 )-plane, 
obtained by Method A. The parameter values are set as tan/3 = 3, A = 0.8, s(0) = 500 
GeV, niA = 1830 GeV, and the temperature is set as T = 100 GeV, which is actually the 
critical temperature T c . Notice two distinct minima of (V(vi, v 2 , T)) on the {v i, w 2 )-plane: 
(0, 0) where the phase of the state is symmetric, and (275, 640) GeV, where the phase 
of the state is broken. The electroweak phase transition may take place from (0, 0) to 
(275,640) GeV on the {v i, t>2)-plane, which is evidently discontinuous and therefore it is 
first order. The distance between the two minima is v c = 696 GeV, indicating that the 
strength of the first-order phase transition is strong (v c /T c > 1). The masses of the three 
scalar Higgs bosons are obtained as ms 1 = 56 GeV, ms 2 = 807 GeV, and rns 3 = 1827 
GeV. 
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FIG. 3: The plot of the equipotential contours of (V(vi,V2, T)) on (t>i, t> 2 )-plane, obtained 
by Method B. The parameter values are set as tan (3 = 3, A = 0.8, s(0) = 500 GeV, = 
1780 GeV, and T c = 100 GeV. The coordinates of two minima are: (0,0) and (165,440) 
GeV. The distance between the two minima is v c = 470 GeV, thus the electroweak phase 
transition between the two minima is strongly first-order. The masses of the three scalar 
Higgs bosons are obtained as = 82 GeV, mg 2 = 804 GeV, and ms 3 = 1777 GeV. 
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